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Abstract: The entanglement entropy in theories with a Fermi surface is known to 
produce a logarithmic violation of the usual area law behavior. We explore the possibil- 
ity of producing this logarithmic violation holographically by analyzing the IR regions 
of the bulk geometries dual to such theories. The geometry of Ogawa, Takayanagi, 
and Ugajin is explored and shown to have a null curvature singularity for all values 
of parameters, except for dynamical critical exponent 3/2 in four dimensions. The 
results are extended to general hyperscaling violation exponent. We explore strings 
propagating through the singularity and show that they become infinitely excited, sug- 
gesting the singularity is not resolved by stringy effects and may become a full-fledged 
"stringularity." An Einstein-Maxwell-dilaton embedding of the nonsingular geometry 
is exhibited where the dilaton asymptotes to a constant in the IR. The unique nonsin- 
gular geometry in any given number of dimensions is proposed as a model to study the 
T = limit of a theory with a Fermi surface. 
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1 Introduction 

Holography is a powerful tool that has recently been employed to study quantum field 
theories that admit a Fermi surface. The renormalization group flow of the quantum 
field theory is geometrized in these models as a flow from large distances to short 
distances in the bulk. 

Another tool to study quantum field theories, which has recently been given a 
holographic interpretation that we will explore in this paper, is the geometric entropy. 
The geometric entropy is defined as the entanglement entropy between two regions A 
and B, where A U B is the configuration space of the entire system. The Hilbert space 
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then admits a decomposition along this separation, H = Ha QQHb- The entanglement 
entropy is computed as 

S A = -Tr A p A \og p A , (1.1) 

where 

p A = Tr B p (1.2) 

with p the density matrix of the whole system. In general the geometric entropy is a 
time- dependent quantity, but in this paper we will focus solely on time-independent 
cases. It was shown that for an arbitrary spatial division of a local quantum field theory, 
the geometric entropy (henceforth referred to as the entanglement entropy) between 
A and B for the ground state density matrix scales as the area of the dividing region 
[1,2]: 

Area(dA) 

oa ~ V subleadmg terms (1.3) 

for ad+l dimensional quantum field theory, e is the UV cutoff of the theory and OA 
is the dividing manifold with Area(cM) ~ L d ~ x for some characteristic length scale L. 
This behavior degenerates in the case d = 1 to a logarithmic scaling: Sa ~ log(£/ e) 
where L is the length of the region A. 

This area law is to be contrasted with the extensive behavior of thermodynamic 
entropy and invites speculation on the connection to black holes. It is intuitively 
understood by observing that for local theories in their ground state, the interactions 
occur most strongly across the boundary. There is a violation of this area law when the 
quantum field theory admits a Fermi surface; in that case, the entanglement entropy 
scales as [3, 4] 

L\ d -\ L 



Sa ~ ^— J log — h subleading terms, (1.4) 

where L is the characteristic length scale of dA. This logarithmic behavior is understood 
as coming from the 2D CFT that approximates the excitations normal to a Fermi 
surface. 

Entanglement entropy has recently been given a holographic interpretation in the 
AdSd +2 /CFT d+1 setup [5]. We imagine a time- independent bulk geometry that is dual 
to some quantum field theory on the boundary. We begin by taking a time slice of the 
boundary field theory. In the time-independent case this corresponds to a time slice of 
the bulk; we work completely within this slice. We split the spatial manifold on the 
boundary of this slice into two pieces A and B with a hypersurface OA. To compute 
the entanglement entropy between these two regions, we first find the minimal surface 
on this time slice in the bulk that is anchored at the boundary on OA. We will follow 
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the convention and call this minimal surface 7^4. Then the entanglement entropy is 
given by 

Sa = (1.5) 

where Gn refers to the bulk Newton's constant. This proposal trivially reproduces the 
"area law," which we show explicitly in Appendix A (for a review, see [6]). 

In this paper, we wish to study the cases in which the holographic entanglement 
entropy formula can reproduce the logarithmic scaling characteristic of the existence 
of a Fermi surface. We will analyze the geometry of [7] (see also [8]), which we call 
the OTU geometry after the authors' surnames, and other cases below. We will show 
that the geometry generically has a null curvature singularity, except for a specific set 
of parameters. The geometry can be modeled by a singular plane wave, and we will 
show that strings in this background become infinitely excited at the singularity, sug- 
gesting that the geometry is singular in the sense of string theory as well. We will find 
that geometries that admit an IR region with a time slice equal to EAdS2 x IR^ 1 are 
the simplest geometries that can reproduce the correct scaling, although their bound- 
ary interpretation is unclear. Geometries such as the electron star and the extremal 
Reissner-Nordstrom black hole will not produce the violation. 

Those who wish to bypass the operatic crescendo of this paper may go directly to 
Sections 3 and 4, which concern the singularities and gravitational embedding of the 
OTU geometry and form the primary results. Some preliminary applications of the 
holographic entanglement entropy formula can be found in Appendix A. 



2 Area Law Violation 

In this section we analyze what criteria are sufficient to violate the area-law scaling 
and produce the logarithmic correction characteristic of Fermi surfaces. We first note 
a seeming problem, which is that any asymptotically AdS^+2 geometry will always 
produce the area-law scaling L d ~ l / e , since this contribution comes from the piece 
of the minimal surface out near the boundary of the spacetime. However, it seems 
clear that the violation should occur in terms of an effective length and an effective 
cutoff (which we call L e $ and e e ff) coming from the piece of the minimal surface probing 
deep into the bulk. 1 This is because in holographic descriptions of Fermi surfaces, we 
expect to begin with a CFT d+1 in the UV of our boundary theory and flow to a theory 

-'^As in [7], we expect e e g ~ kj, 1 (this parameter will sometimes be set to 1) for the Fermi momentum 
hp, so statements concerning the IR geometry arc meant to imply the limit z 3> e ff- 
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with a Fermi surface in the IR. The IR of the boundary theory corresponds to regions 
deep in the bulk geometry, so the logarithmic violation should be produced there. Of 
course, there are subtleties in this interpretation, but we will adopt it as our operating 
assumption to see what answers we get. 

Instead of working in the full geometry and considering large minimal surfaces 
that probe deep into the geometry, we can analyze the IR region of the geometry 
alternatively by pushing the cutoff deep into the bulk and excising the UV region. 
Then, the geometry on our boundary theory slice will be different and we can hope 
to get a different type of scaling as we approach this new "UV" boundary; this will 
be made precise in an example below. We will find that a geometry that becomes 
AdS3 x IR^" 1 deep in the bulk (or, more generally, its time slice becomes EAdS2 x IR d_1 ) 
will satisfy the logarithmic violation with respect to L e g and e e ff. 2 This is simply because 
the AdSs piece will give the logarithm we are looking for, while the piece will give 
the "extensive" L^' 1 / e^ 1 piece. They will combine to give (Lfg 1 j efg 1 ) log (L eS / e cff ), 
which is the result we hope for in an asymptotically AdSd +2 geometry. We will also see 
that the electron star and extremal Reissner-Nordstrom black hole geometries do not 
give the correct scaling. Adopting a certain metric ansatz, the authors of [7] produce 
general geometries which contain the logarithmic violation; we will present the geometry 
and analyze its singularity nature in Section 3. 

Viewed in this way, working in the IR geometry is simply a shorthand to extract the 
IR pieces of the full minimal surface. This approach misses important UV-IR mixing 
terms since the UV part of the geometry is excised. It is possible that logarithmic 
violations can come from these intermediate regions. 

When working in the IR geometries, we will consider them geometries in their own 
right and apply the holographic entanglement entropy prescription from their "bound- 
aries." When we use this language it is simply shorthand for the procedure involving 
pushing the boundary slice deep into the IR region. This distinction is important be- 
cause the IR boundaries we consider are not meant to be considered as holographic 
screens (or as conformal boundaries), and as stated above serve only as a tool to ex- 
tract the IR contribution to the entanglement entropy. 

2.1 AdS 3 x R^ 1 

In this section we will analyze 3 AdS3 x R d_1 (again, any geometry with a time 

2 e ff refers to the new location of the boundary slice, while L e g refers to the characteristic size of 
dA on this slice. See Appendix B for an example of how we expect e e ff and L e g to arise from the full 
geometry. 

3 The entangling properties of this spacetime were also studied in [9] 
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slice of EAdS2 x IR d_1 will do the trick, but we stick to this geometry for concreteness) , 
which we hope to appear in the IR of some asymptotically AdS^+2 geometry We write 
this IR metric as 

ds m = \(~dt 2 + dz 2 + dy 2 ) + \-dx\ (2.1) 

for i — 1, 2, . . . d— 1 and i?Ads = 1- We follow the warping to the regulated "boundary" 
at z — e c g. Since the IR metric is a simple product manifold, the minimal surfaces will 
be product manifolds as well. We already know the AdS3 answer for a minimal surface 
(in this case just a spacelike geodesic); for a length L c s on the boundary it is given by 

(z,y) = ^(sins,coss), (2.2) 

which has length 

Length ~ log — — . (2.3) 

e c fr 

For the full boundary geometry including the ]R d_1 piece, we consider a d— 1 dimensional 
box of side length L e g separating a time slice of the boundary into regions A in the 
interior of the box and B representing the exterior. Restricting to the part of the 
geometry, this becomes a d — 2 dimensional box of side length L e q. The minimal surface 
anchored at the boundary on this d — 2 dimensional box is exceedingly simple to find; 
you simply fill in the interior of the box! Since no warping appears, the surface remains 
on the boundary slice and does not reach into the bulk. The area of this minimal 
surface scales as 

Area-^, (2.4) 

6 eff 

which is simply the extensive scaling observed in Appendix A. 3. Considering the prod- 
uct of these two minimal surfaces, we get the appropriate scaling 



Area- ^ log ^. (2.5) 
e eff e efr 



Notice that all we required in this analysis is that a time slice of the geometry be 
EAdS2 x ]R d_1 . For example, modifying to Lifshitz scaling will not change these results, 
as long as it is a Lif 3 piece analogous to our AdS3 piece, which will keep the same induced 
geometry on a time slice. 

The physical interpretation of the boundary theory is less clear. The AdS3 factor 
suggests that the specific heat (and entropy density) depend linearly on temperature, 
just as a conventional Fermi liquid. A troubling aspect of this geometry is that it is 
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usually produced by turning on magnetic fields, in which case the low energy theory be- 
comes effectively 1+1-dimensional (describing the physics of the Lowest Landau Level), 
and the logarithm we are picking up is simply that of the 1+1-dimensional theory. Such 
a model is given by the magnetic brane solution of [10, 11]. At this level, this geometry 
is presented simply as a proof of principle for the existence of the logarithmic violation, 
before moving on to the more general geometries of [7]. 

2.2 Other Fermi Surface Geometries 

In this section we analyze some of the other geometries that are purported to have 
Fermi surface duals. One model is the extremal Reissner-Nordstrom black hole [12-15] 

R 2 f dz 2 \ 

ds 2 = tL (~f(z)dt 2 + + dx 2 j , (2.6) 

where 

/w=i - 4 (£) 3+3 (£f (27) 

is the emblackening factor. Expanding this geometry in the near horizon limit (around 
z = zq) and rescaling coordinates one gets 

R 2 

ds 2 = -^(-dt 2 + dp 2 ) + dx 2 , (2.8) 

which is simply AdS2 x W 1 . The induced metric on a constant time slice can be written 

as 

R 2 

ds 2 = -^dp 2 + dx 2 = dp 2 + dx 2 , (2.9) 

which is simply R +1 . This will produce an extensive entanglement entropy contribution 
Lf s / ef ff as argued in Appendix A. 3. The fact that the minimal surface doesn't reach 
into the deep IR (the "boundary" is placed at p — e e fr and any minimal surface in the 
bulk anchored at the hypersurface dA on the boundary actually lives completely on 
the boundary slice) mirrors the fact that minimal surfaces in the full geometry tend 
to hug the horizon. This gives the extensive behavior shown in [16, 17] by numerical 
integration in the full geometry. 

The electron star is similar but in the near-horizon limit has a Lifshitz form defined 
by a different scaling for the coefficient gu [18] 

ds " = R2 (-^dt 2 + 9 -^dr 2 + ^dx 2 ) , (2.10) 
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where g^ is a constant that depends on various parameters in the construction and 
z is the dynamical critical exponent. Since we work on constant time slices, the new 
coefficient gu is irrelevant, and we see that a constant time slice in the near horizon 
geometry is simply EAdSd+i- This will give the area- law scaling and again no logarithm. 

2.3 The OTU Geometry 

A geometry with the near-horizon limit AdS3 xlR^ 1 is a simple situation in which a 
logarithmic violation can occur. However, it is possible to have geometries that do not 
factorize in the IR but produce the correct entanglement entropy scaling for a Fermi 
surface. 4 We shall present the asymptotically AdS4 geometry in [7] and analyze its 
singularity nature in Section 3. 

One begins by picking a metric ansatz 

B 2 

z 2 



ds 2 = — (-f(z)dt 2 + g(z)dz 2 + dx 2 + dy 2 ) (2.11) 



in four dimensions and requiring the logarithmic behavior with an effective cutoff. Ac- 
cording to [7], this fixes the function g(z) to z 2 / e 2 s deep in the bulk. Since we consider 
time- independent geometries, f(z) is not fixed by the behavior of the entanglement en- 
tropy and is instead partially fixed by imposing the null energy condition in the bulk. 
We ignore f(z) for now and focus on a given time slice. 

As usual, we go to the "boundary" of the IR geometry (which we identify as z — > 0) 
and pick a region A defined by the infinite strip: 

A := x E (-L x /2, L x /2), y E (~L y /2, L y /2) (2.12) 

with L y ^> L x . As usual, these are effective lengths L e g seen in the IR when embedded 
in the full geometry, but we leave off the subscripts to keep the equations less cluttered. 
The area of the minimal surface is given by 



1/1 / z' N " 



Area( 7A ) = J dxdy-d — + |_J . (2.13) 

Translation invariance in y due to the belt geometry allowed us to parametrize the 
surface as z(x,y) = z(x), with z' = dz/dx. We can then perform the ^/-integral to get 



*W2 i i / z i \ 2 
Area( 7A ) = L y I dx-J- + — . (2.14) 

-L x /2 Z V Z \ e eff/ 



4 This analysis was partially carried out during the course of this paper, but a fuller treatment was 
given in [7]. 
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Considering the integrand as a Lagrangian £, then 

*-£-<-'-^(H£fT m -* (215) 

is a conserved quantity, where we imagine x as a time variable. The minimal surface 
will be a symmetric function of x, which gives z'(x = 0) = for the turning point. We 
use this conserved quantity to write 

z 3 dz , f a a 2 dz 

Area = 2L V t^=^ ( 2 -16) 



e cS e eS a/o 4 - z 4 v J eeS e eS z^ a 4 - z 4 

where we cut off the integrals to extract the divergences. These integrals can be com- 
puted analytically and give 



-eff 



a 2 + v / c 



err 



Area = L y ^ '- . (2.17) 

e c fr e e fr 

This is one of those rare and delightful cases where expanding for small e e g is easier 
without the aid of M athematica: 



L x ~ — Area ~ L y (2.18) 

Ceff e eff 



a 2 



Area ~ hi l og ( 2 ^-) ~ h. l og fe^) , (2.19) 

where we have ignored terms subleading in e e g. A similar analysis can be applied to 
higher dimensional cases. We see that our IR prescription produces the correct scaling 
in the geometry identified in [7]. 



3 Singularities of the OTU geometry 

We wish to analyze the singularity nature of the OTU geometry. In the original 
paper [7], the authors find an effective gravity embedding of the geometry with fields 
which become singular at z = oo. We will show that the geometry generically has a 
null curvature singularity at z = oo except for a specific set of parameter values. The 
singularity exists even though all curvature invariants remain finite. In fact, the sin- 
gularity nature is analogous to that of the Lifshitz geometry, in a sense which we will 
make precise below. The method of analysis is inspired by [19], and we will compare 
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our results to those of the Lifshitz geometry, which is included as a special case of our 
geometry. 



3.1 Diverging Tidal Forces 

We begin by considering the metric 

.2 1 



ds 2 = — (-f(z)dt 2 + g(z)dz 2 + dx 2 + dy 2 ), (3.1) 

z 2 

with f(z) ~ z~ 2m and g(z) ~ z 2p . The radius of curvature is set to 1. We will restrain 
from specifying to p = 1 and m > p, the conditions required in [7], until the end. We 
also note the alternative form of writing the metric which makes the dynamical critical 
exponent 7 and hyperscaling violation exponent 9 manifest, for later comparison [8, 20]: 

* 2 = p + + **) ■ < 3 - 2 > 

Comparing to our metric shows that d = 2, m = 2(7 — l)/(2 — 6 1 ) and p = 0/(2 — 9). 

We want to consider the tidal forces felt by an infalling observer, so we consider 
an orthonormal frame that is parallel propagated along a radial infalling geodesic. The 
four- velocity is given as £/ M = (t, z, 0, 0), where dots denote derivatives with respect to 
proper time. We have as an integral of the motion the conserved energy E = —K^U^ = 
—Qtt%- an d the normalization of the four- velocity gives 

z 2 ( E 2 z 2 ' 



z 2 = ^-(-l + —^). (3.3) 
9(z) V f{z)J 

We write the vielbein parallel propagated along a radially infalling geodesies as 

(eof = Ez 2m+2 {d t Y + zW-l + E 2 z 2 + 2m {d z y (3.4) 

(e x Y = -z m+ W-l + E 2 z 2+2m {d t Y - Ez 2+m - p {d z y (3.5) 

(e t y = z(d x .y. (3.6) 



Components of the Riemann tensor in this frame will give the physical tidal forces felt 
by an infalling observer. We use the constructed vielbein to transform from the "static" 
frame to the infalling frame: 

Robed = R»u P v(e a y(e b y(e c y(e d y, (3.7) 
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and we find the following diverging components in the new frame 




2p (l+p + E 2 (m-p)z 2+2m ) 
2p (-l-m + E 2 (m-p)z 2+2m ) 



(3.8) 
(3.9) 
(3.10) 



E{m - p)z 1+m - 2p y / (-l + E 2 z 2 + 2m ). 



Restricting to m > p > 0, the tidal forces diverge as (m — p)z 



-2 P +2_ Notice that 



for p = we simply get the result in [19] which is appropriate for Lifshitz scaling. 
m = p = is the nonsingular result of pure AdS. However, in this more general 
geometry we see that there seem to be other nonsingular results, namely for m = p 
and m < p — 1. Imposing the null energy condition in the bulk gives the condition 
m > p, which leaves only m = p as a nonsingular set of values for the parameters. 
This constraint translates, in general dimension, to 7 = 1 + 9/d for the condensed 
matter exponents. In general dimension this also saturates the null energy condition 
in the bulk. To reproduce the logarithmic violation, we need to fix p = 9 = 1 in four 
dimensions and so we are left with the nonsingular metric 5 



The hyperscaling violation exponent for this geometry is 6 = 1 while the dynamical 
critical exponent is 7 = 3/2. Surprisingly, these values occur precisely for gauge theories 
of non- Fermi liquid states (in a three- loop analysis) [21, 22]. This metric has a Poincare- 
type horizon at z = 00, to which timelike geodesies can propagate in finite proper time 
(one can see this by integrating (3.3); the singularities in the singular geometries are 
also reached in finite proper time). But the metric should admit an analytic extension of 
the coordinates, just like the extension from Poincare coordinates to global coordinates 
in AdS. A quick computation of the surface gravity shows that z = 00 is a degenerate 
Killing horizon for the Killing vector dt- Some intuition for the special status of this 
metric can be gained by making the coordinate transformation z = 1/r, after which 

5 The translation invariance of this metric in the coordinates x\ and X2 suggests that geodesies 
with momentum in these directions, not considered above, should not produce divergent tidal forces 
(although this has only been checked explicitly for i?0303 with momentum in the 2-direction) . Further- 
more, curvature invariants will all remain hnite at z = 00, since in the "static" frame all components 
of the Riemann tensor are finite there. The singular cases are similar to the Lifshitz spacetime, in 
the sense that curvature invariants are all finite while there exists a null curvature singularity. The 
singularity at z = is of no importance because this geometry is only valid for z 3> kp 1 = 1; it is 
asymptotically AdS, which is well-behaved at z = 0. 
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the metric takes the pleasant form 



ds 2 = -F(r)dt 2 + -^-y + r 2 dx 2 , F(r) = r 4 (3.12) 

The horizon is now located at r = and we see that a continuation to negative r 
preserves the signature of the metric. 

The singularity for m > p > is a null curvature singularity, as surfaces of constant 
z become null as z — > oo. This can be seen by observing radial null geodesies: dt = 
±z 2m+2p dz. Since our singularity is at z = oo, it takes an infinite amount of coordinate 
time to reach the singularity, implying that the singularity is null (recall we work in 
the regime m > p > 0). 

In higher dimensions, as long as rotational and translational symmetry is kept in 
all dimensions but z and t, the results are similar. The value of p will increase, but the 
null energy condition will still give m > p. 

3.2 Plane Wave Approximation 

In this section we show that, as in the Lifshitz geometry, the OTU geometry can 
be represented near the singularity as a plane wave. To do this analysis we work with 
the geometry in the form 

ds 2 = {-r 2m+2 dt 2 + ^- 2 + r 2 {dx\ + da*)) , (3.13) 

where we have simply made the redefinition r — 1/z. Switching to a tortoise coordinate 
via dr*/dr = y/—g rr / gtt =^ = — r n ^_ p P +1 and subsequently lightcone coordinates 
u — t — r*, v = t + gives the metric in the form 



2m+2 2 

ds 2 =(^r) ° (-^ 2 + ^)+(^) C (dxl + dx 2 ) 

m+l I 

dudv + I — r- I (dx\ + dxl), 



c 2 (u — v) 2 J \c 2 (u — v) 



(3.14) 



where we have defined c = m + p + 1. Since r ~ (r*) -1 ^ ~ {^z^) 1 ^ \ we have that 
small r corresponds to u ^> v if we stick to c > 0. The viable set of parameters for 
logarithmic violation of the entanglement entropy is in this range (m > p = 1 in four 
dimensions). Thus, near the singularity we can approximate the metric as 
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m+1 

4 



If we make the redefinition 



c 2 u 2 



du = dU, (3.16) 



then for the coefficient of dx 2 + dy 2 we get 

l/c 



(i) ^ = (3.17) 



which gives a plane wave metric 



z 



2 

Notice in the case p = 0, 2m + 2 = 27 for critical exponent 7, we have precisely 
the Lifshitz metric and the answer we get is precisely the Lifshitz answer, with the 
coefficient of the second term being given as (jU/2) 2 ^. In general, we can write our 
oefficient in the suggestive form (kU/2) 2 ^ k for k = m + 1 — p\ this is precisely the same 
form as the Lifshitz answer! 

To bring the plane wave metric into the standard Brinkmann form, we define 



X X i /kU\ l ^ k 

k = m + l-p, V = V-^—, Xi = X A — \ (3.19) 



and get 



ds 2 « -dUdV + dXidX 1 + W^X^dU 2 , W{U) = \— \. (3.20) 

This is a good approximation to the geometry near z = 00 as long as k ^ 1, since the di- 
verging tidal forces of this geometry and the OTU geometry behave in the same way. If 
k — 1, however, the OTU geometry is some nonsingular, asymptotically AdS geometry 
while the plane wave approximation gives Minkowski space. Thus, the approximation 
breaks down in this there are no diverging tidal forces which dominate the 

behavior and match up between the two geometries. The metric is appropriate in any 
number of dimensions, where p will change to give the logarithmic violation. In all 
cases, however, p = m ==>■ k — 1 will be the one nonsingular case which also satisfies 
the null energy condition. 
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The calculation of infalling geodesies is identical to that of the OTU geometry. The 
four- velocity is given as £7 M = (If, V, 0,0), where dots denote derivatives with respect 
to proper time. We have as an integral of the motion the energy E = —K^U^ 1 — U/2. 
Normalization of the four-velocity gives 

l + ^-W(U)X,X' 

2E K ' 

We pick the parallel propagated frame 

(e»r=2E { a u r + 1 ±^MMi {9vr (3 . 22) 

( tf = 2gW+ ±iM! W (3.23) 
(e i r=(d x >r. (3.24) 

We find diverging components of the Riemann tensor 

Rom = Run = Rom = AE2 ^u2 l) = ( m " P)z 2(m+1 - p \ (3.25) 
which diverges precisely like the original metric as z — > oo. 



3.3 A Stringularity 

Studying test strings in the plane wave background (3.20) is a solved problem; 
the strings propagating through the singularity become infinitely excited, suggesting 
that the geometry is singular in the string theoretic sense, i.e. "stringular." However, 
without a string theory embedding this should be viewed merely as a hint; as is well 
known by now there are other effects, e.g. the dynamics of D-branes, which can resolve 
singularities. 

The mathematics of these test strings is identical to [19], which in turn is essentially 
identical to [23]. We will calculate the expectation value of the mass squared and mode 
number operators to see that they diverge as long as m ^ p. This more general case for 
the value of p is relevant for higher dimensions, as imposing the null energy condition 
in the bulk always gives m > p. 

We begin with the usual Polyakov action describing strings propagating in a general 
curved background 

S = --^ / drda^^g^daX^d^. (3.26) 



13 



We can fix some of the gauge redundancy by working in conformal gauge = e^ a ' T ^r] a b. 
The residual freedom (changing a and r by solutions to the wave equation) is fixed by 
imposing light cone gauge, U = a'pr. Fourier expanding the embedding coordinates as 



n=+oo 



xw)= E x n^v na ( 3 - 27 ) 

n=— oo 

allows us to write the equations of motion in this gauge as 
X* n + [n 2 - a' 2 p 2 W{a'pr)} X l n = 

a'pV = {X' 1 ) 2 + (X H ) 2 + a r2 p 2 W{a' F r)XX i (3.28) 
a'pV = 2X i X'\ 

where we see that the equation of motion for the embedding coordinates is a Schrodinger 
equation. 6 From this equation we see that the mode functions depend on the combi- 
nation nr: 

Xi(r) = Xi(nr). (3.29) 

We consider a plane wave with flat spacetime before and after it. In other words, we 
take 

, x fo if \U\ > T . 



1 -=k ( A L.) if 1 7"/ 1 < T 



fc 2 



for some large T, where is negative for the singular cases of interest (m > p), giving 
an attractive potential in our Schrodinger equation. For the regions \U\ > T we can 
do the standard flat space expansions of the embedding coordinates: 

X\a,r) = q i < + 2a'f <T + ^e»i; < (r) r < -r (3.31) 

X\a, r) = ql + 2a'p i > r + ]T e^X^r) r > -r (3.32) 

K>,<(r) = ^«>,<e- mr - 5 l _ n>i< e^) n ± (3.33) 

6 This Schrodinger equation has potential a' 2 p 2 W(a'pT) = ^ra near the singularity. Notice that, 
contrary to the statements in [23] which ignore the strength of the potential, for the cases of interest 
k > 1 and the strength of this inverse square potential is bounded between —1/4 and 0, allowing for 
transitions in the quantum mechanics problem; thus, the string can make it through the singularity 
at t = and our matching is justified. 
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for t = T/a'p and mode operators a l n> <1 a l n> < which satisfy the canonical commuta- 
tion relations. The slightly unconventional normalization has been chosen to conform 
to the past literature on the subject. The Bogoliubov transformations connecting these 
two expansions are written generally as 

<> = K<< + K4 f < (3-34) 

h { - 4* n 1 -I- W n*t 
u n> — J-i n u n< t u n u, n< . 

The Schrodinger equation for X l n< can be solved exactly by Bessel functions and 
expanded as r — > CT to get 

X l n< (r) = D* Inrp-")/ 2 + F^nrf^' 2 for 1/ = ^1 - 4C(1 - k)/k. (3.35) 

The problem is essentially a scattering problem off of the given potential. Accord- 
ingly, we pick as an initial state 

f n {r) = for r < -r (3.36) 

and get a solution of the equation of motion for the embedding coordinate as 



2 12 

fn(r) = ^ nr ^~ ( e mr / dr'e- inT 'r n (T')W(a'pr') - e~ mT \ dr'e™' ' f n (r')W (a' pr' 



(3.37) 

One can verify this is a solution by plugging into (3.28). From the equations of motion 
for our embedding coordinates (3.31), (3.32) and the Bogoliubov transformation (3.34) 
we get 

K = S- dTeinT ^ T ) w ^'P T )- ( 3 - 38 ) 

The mass and number operators in the U > T region are given by 



Ml 



\, ( a n> a n> + + m l> (3-39) 

n=l 

oo . 

- (4><> + 4>4>) > (3-40) 



oo 

iV> = 

n=l 

where Wq = — (£) — 2)/(12a') is the tachyonic mass in this normalization. 

Using the Bogoliubov transformations (3.34), we can find the expectation values 
of the mass and number operators in the ingoing ground state as 

x <l </ n=i i=i 

' ^ ' < ' n=l i=l 
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for D spacetime dimensions. Making the change of variables y = nr in (3.38) and using 
(3.29) and the fact that T is large gives 



2%k 2 J_ nT0 y 2 2k 2 7_ nro y 2 

This integral is finite and independent of n for k ^ 1 [19]. Thus, each mode is excited 
equally and mass and number operators diverge, leaving the door open for a full-fledged 
stringularity. We are only saved when k — 1, which corresponds to m = p. This is 
precisely the case where there is no general relativistic singularity. We stress again that 
this calculation is intended to be suggestive rather than a definitive remark about the 
behavior of this geometry in string theory. 



4 Effective Gravity Embedding of Nonsingular Geometry 

The OTU geometry for the singular m > p = 1 case is derived from Einstein- 
Maxwell-dilaton theory in [7]. We now demonstrate a solution in the full geometry 
for the nonsingular case m = p = 1. The geometry we want to produce is ds 2 = 
ji (-f(z)dt 2 + g(z)dz 2 + dx 2 + dy 2 ) with f(z) = 1/(1 + z 2 /k) and g(z) = 1 + z 2 /z 2 F . 
This has the required nonsingular IR behavior which produces the logarithmic violation 
of entanglement entropy, in addition to being asymptotically AdS±. It satisfies the null 
energy condition everywhere in the bulk when k < Zp, which is necessary for the reality 
of the scalar field solution given by (4.4). 

We begin with the action 



S 2k 



(R - 2A) - Z{4>)F^F^ - \d^d^ - V(<P) 



(4.1) 



With an ansatz = <J)(z) and A t = a(z), we get a'(z) = -^^y/ f(z)g(z) for an integra- 
tion constant A. All equations of motion can be solved by 
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v(4>) 

z 2 g{z)f{zf - Af{z)\Qg{z) - 6g(z) 2 + zg'(z)) + zf(z)(zf(z)g'(z) + g(z)(8f'(z) - 2zf"(z))) 



Z(cf>) 



±f{z) 2 g{zY 

U 2 z 3 f(z) 2 g(z) 2 



zg{z)f'{z) 2 + f(z)(zf'(z)g'(z) + g(z)(4f'(z) - 2zf"{z))) : 



zf{z)g(z) 

We place this preliminary form to highlight the ease with which metrics can be produced 
with this action and ansatz. For the specific f(z) and g(z) proposed above, we find 

Mz) = AZF ^ (4.3) 

2z{k + z 2 fl 2 ^z 2 + z 2 F 

<p(z) = 2yjl - k/z 2 F F axctm(z/Vk) 1 - k/z 2 F , (4.4) 

where F[a|m] is the elliptic integral of the first kind. Constants q have been defined to 
condense the expression; they are just combinations of Zp and k. The dilaton is plotted 
in Figure 4. We can invert this to get z in terms of (f) and get expressions for V(4>) and 



_ sc[0/(2 v /cl);ci ] 2 (c 2 + c 3 sc[0/(2^cl); d] 2 + c 4 sc[0/(2 v /cl); Ci] 4 + c 5 sc[</>/ '(2^/cT); c x ] 



1 + sc[</>/(2^c7); Cl ] 2 ) 2 (4 + fec[0/(2^7); c x ] 



\2\2 



2fc 2 sc[0/(2y^); c] 2 (l + sc[0/(2^); c] 2 ) 2 (z 2 F + fcsc[0/(2^); c] 2 ) 2 
{<P) ~ z 2 F (z 2 F + 2{k + 2z\ )sc[0/(2V5); c] 2 + 5tec[0/(2VH); c] 4 ) ' 1 ' ] 

where sc[a; b] is the Jacobi elliptic function, which appears due to inverting the elliptic 
integral in (4.4). These are not illuminating forms of the functions; their more tangi- 
ble plots can be found in Figures 2 and 3. We see that the dilaton asymptotes to a 
constant, whereas in the singular cases where m > 1 it diverges logarithmically. As 
one can guess, this constancy of the dilaton in the IR is robust to modifications of the 
metric which keep the asymptotic behaviors the same. 
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3.0 



4(4+2(fe+24)sc[0/(2^);c] 2 +5fcsc[0/(2 v / S);c] 4 ) 



with k = 1 and zp = 2. 



The divergence is simply indicating that the gauge coupling, which scales as Z(4>) , is 
vanishing in the IR. 

a(z) 
6 



3 - 



1 




2z(k+z 2 ) s / 2 y/z 2 +zj 
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5 Conclusions 

In this paper, we analyzed situations in which the logarithmic violation of the 
entanglement entropy, indicative of a theory with a Fermi surface, is expected to occur 
holographically. This analysis uses weakly coupled gravity in the bulk, which means we 
are assuming that the logarithmic violation persists to strong coupling. The AdSsxIR ' -2 
geometry was presented as a simple case which produces the violation, though the dual 
interpretation is unclear. 

The geometry of [7] was presented as a more general form which can produce the 
logarithmic violation. It has gained support from analyzing the scaling of the thermal 
entropy density and it has been argued that the holographic entanglement entropy 
formula only detects "hidden" Fermi surfaces [8] . Further holographic explorations are 
performed in [20]. This class of geometries generically has a null curvature singularity 
at z = oo which is reached in finite time, except for when m = p. In this case, the 
surface z = oo is like the Poincare horizon of AdS, and the metric can be analytically 
extended beyond. In four dimensions, this metric takes the form 



We can gain further intuition for this candidate by writing down the general metric 
with dynamical critical exponent 7 and hyperscaling violation exponent 9 [8, 20] 



The logarithmic violation of entanglement entropy occurs for 9 = d — 1 which equals 1 
in our case of 4 bulk dimensions. We then see that 7 = 3/2, which saturates the bound 
from the null energy condition. It is interesting to note that these specific values of the 
exponents are what arise in a three-loop analysis of gauge theories of non-Fermi liquid 
states [21, 22]. More generally, the nonsingular cases are given by 7 = 1 + 9/d, which 
again saturates the null energy condition. 

The singular cases are only an obstruction to studying the far IR of the theories, 
as there exist black hole solutions which mask the singularities; in this sense, they are 
"good" by the classification of Gubser [24] and the limit of T — > can be taken, which 
corresponds to the finite temperature IR cutoff being removed in the dual theory. How- 
ever, to study the T = point, the solution with in = p should be employed to avoid a 
naked singularity. An embedding of this geometry into Einstein-Maxwell-dilaton the- 
ory has been provided in Section 4. The action considered does not permit a consistent 




(5.1) 




(5.2) 
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IR truncation (because the dilaton asymptotes to a constant), but if a different action 
with such a consistent IR truncation is found, then black hole solutions with specific 
heat C ~ T 2 / 3 are guaranteed to exist [7, 8, 20, 25, 26], allowing you to study the 
system consistently both at zero and finite temperature. 
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A Holographic Entanglement Entropy Computations 

In this appendix we apply the proposal to compute the entanglement entropy in a CFT 
holographically. We show how the area-law scaling comes naturally from the asymptot- 
ically AdS geometry. We begin with the logarithmic scaling in three bulk dimensions 
and move to higher dimensions. We end with a consideration of the extensive scaling 
in flat space. 



A.l Entanglement Entropy in 2D CFT 

We will work in Poincare coordinates 

R 2 

ds 2 = —(-dt 2 + dz 2 + dx 2 ). (A.l) 

z l 

The boundary theory is 1+1 dimensional, so the region A is an interval of length L 
at a fixed time slice t = to. It is divided from region B by the points x = —L/2 and 
x = L/2. The minimal surface in this dimensionality reduces to a spacelike geodesic. 
To get finite answers we put the boundary theory at the slice z — e. The geodesic ja 
can be parametrized by the non-affine parameter s as 

(z, x) = — (sin s, cos s). (A. 2) 

Computing the length of this geodesic, we get 

Length^) = / J g .. ds= — = 2i?log-. (A.3) 

J eq » us as J eo sin s e 
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where e = 2 e / L and is the induced metric on a constant time slice. This gives for 
the entanglement entropy 

This reproduces the logarithmic scaling we expected. 

We see in this computation the famous UV-IR connection; an IR cutoff in the bulk 
(picking a slice a finite but large distance away, z — e) corresponds to a UV cutoff in 
the boundary theory. 



A. 2 Entanglement Entropy in Higher Dimensions 

We work again in the Poincare coordinates 

ds 2 = ^j(-dt 2 + dz 2 + dx 2 ) (A.5) 

for i — 1,2, . . . ,d — 1. In higher dimensions there are more shapes we can pick for the 
region A on the boundary z — e. The computation of the minimal surface simplifies 
when a region is picked that maintains as much of the symmetry as possible. We will 
pick the simplest such surface, the "belt" geometry, which has x\ G (— L/2, L/2), Xi G 
(— L±/2, L±/2) for % — 2, 3, . . . , d — 1 with L± 3> L. In this case, we minimize the area 
functional 



1/2 y/1 + (dz/ds) 
L/2 



2 



i^Lf 1 / ^ — J -ds (A.6) 



and get 



R d ( L 



2G N (d-l) 
which is the expected area law behavior. 



d-i 

+ subleading terms, (A.7) 



A. 3 Comparison to Flat Space Scaling 

Here we consider the scaling of entanglement entropy in flat space. The minimal sur- 
faces are simple to understand, and one can quickly see that for a given slicing of flat 
space, one gets extensive scaling of the minimal surface, i.e. the minimal surface in 
the bulk anchored at dA on the boundary scales as the volume of OA as opposed to 
its area. This can be simply understood to be coming from the fact that a minimal 
surface connecting a large rf-dimensional circle (the shape is not relevant; we pick a 
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circle for concreteness) in flat space is a d + 1-dimensional ball, which scales exten- 
sively. We contrast this extensive scaling Sa ~ L d / e d in ]R d+1,1 with the AdS^+2 result 
Sa ~ L d ~ l j e d ~ x . This distinction will be important below. 



B Working in the Full Geometry 

A satisfying way to investigate whether a logarithm appears is to work in the full 
geometry and pick large surfaces that probe deep into the bulk. In this case, the 
area of this large surface will come with the usual UV divergent area-law piece (which 
accompanies any asymptotically AdS spacetime) but will have additional pieces that 
correspond to the surface reaching deep into the bulk. It is possible to identify an 
effective cutoff and effective length that appear, e e fr and L e ff> and look for a piece 

Td-l T 

-f=r lo s — ■ (B.i) 

e eS e eS 

We illustrate this behavior with an illuminating toy example constructed in [27]; there, 
the authors consider a bulk geometry consisting of two AdS pieces patched together 
with a sharp domain wall 

ds 2 = e 2A{r \-dt 2 + dx 2 ) + dr 2 } (B.2) 

where 

A(r) = l r/RuV **r>r im (B.3) 
[r/R IR for r < r DW 

This is a discontinuous version of realistic models with a continuous domain wall, like 
the RG flow from M = 4 Super- Yang Mills to an M = 1 fixed point (called the Leigh- 
Strassler fixed point [28]) triggered by a mass term for one of the chiral multiplets. 
Continuing with the toy model in AdSs, we find a minimal surface anchored on the 
boundary on a ball of radius I. The area of the minimal surface is found to be 

- \e 2r ™l R ™R uv lJp - \ + l -Rl v tanh" 1 
+ ^e 2rDw/RlR \/l 2 - lljl 2 - I 2 + R 2 e- 2r ^/Rm 




2 V <•;• V • 'T IR 




^tanh" 1 I J - 11 I +<D(e), (B.4) 

1 1 " -l^ + R^e-^w/Rm ' 
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where the definition 

ll = R 2 uv e- 2r ™' Rvv = l 2 cr + 0{e 2 ) (B.5) 

has been made. The first line in (B.4) would be the answer if we were in pure AdS. 
There are also purely IR terms and IR-UV mixing terms. Identifying 

L 2 cS = l 2 ~ ll + 0(e 2 ), e eS = R IR e- 2r ^ R '» (B.6) 

and expanding the terms in (B.4) proportional to Rj R for e e s / L e s <C 1, the authors 
get an IR contribution 

Hp (ff + l0 S (g)) + % X - 2 M2)) + 0(e cS /L cS ), (B.7) 

which is the IR version of the pure AdS answer. We thus see how an effective length 
and cutoff appear in terms of sensible parameters of the theory and give the expected 
answer in the IR. In the realistic models, the minimal surface continues to break up in a 
sensible way, with the piece on one side of the domain wall giving information relevant 
to IR physics and the piece on the other side of the domain wall giving information 
relevant to UV physics. 

This toy example is an illustration that the minimal surface of the full geometry 
is smarter than it looks. When we excise the UV geometry and work solely in the IR 
geometry in our analysis, we are imagining our answers as being the IR pieces coming 
from a more complete procedure like the one above. Of course, important UV-IR mix- 
ing terms are not picked up by this analysis. 
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